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Abstract 

The hard-gluonic contribution to the first moment of the polarized proton structure function 
Qi{x) is dependent of the factorization convention chosen in defining the quark spin density and the 
hard cross section for photon-gluon scattering. Two extremes of interest, namely gauge-invariant 
and chiral-invariant factorization schemes, are considered. We show that in order to satisfy the 
positivity constraint for sea and gluon polarizations, the polarized valence quark distributions 
should fully account for the observed g\{x) at x > 0.2. This together with the first-moment and 
perturbative QCD constraints puts a pertinent restriction on the shape of Au v (x) and Ad v (x). 
The spin-dependent sea distribution in the gauge-invariant factorization scheme is extracted from 
the data of g\{x). It is shown in the chiral invariant scheme that it is possible to interpret the 
g\(x) data with anomalous gluonic contributions, yet a best least x 2 fit to the data implies a gluon 
spin distribution which violates the positivity condition |AG(x)| < G(x). We then propose a more 
realistic set of parton spin distributions with sea polarization and with a moderate value of AG. 
The polarized parton distributions in this work are presented in the next-to-leading order of QCD 
at the scale Q 2 = lOGeV 2 . Predictions for the polarized structure functions fff (x) of the neutron 
and gf(x) of the deuteron are given. 
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I. Introduction 

In the last few years we have witnessed a remarkable progress in the study of polarized 
hadron structure functions and the related proton spin issue. Experimentally, new measure- 
ments of the longitudinal spin-dependent structure functions on various targets in polarized 
deep inelastic lepton-hadron scattering became available. The polarized structure functions 
<7i(x) of the proton [1,2], g'±(x) of the neutron [3], and gf(x) of the deuteron [4,5] have been 
measured recently. The original EMC experiment on g\{x) [6], which has triggered a great 
deal of interest in the proton spin structure, is confirmed by the new high-statistics exper- 
imental data. Theoretically, a direct first-principles lattice QCD calculation of the proton 
matrix elements of the axial vector current, which is free of the r/' and related problems 
encountered before [7], is also available very recently [8,9]. The calculated quark spin is 
consistent with experiment. It is also evident from the lattice calculation that it is the dis- 
connected diagram, which is presumably dominated by the axial anomaly, that explains why 
the total spin carried by the quarks in a polarized proton is smaller than naively expected. 

In spite of the aforementioned progress, the extraction of spin-dependent parton distri- 
bution functions, especially for sea quarks and gluons, from the measured polarized hadron 
structure functions remains largely ambiguous and controversial. One main issue has to 
do with the debate of whether or not gluons contribute to 1^, the first moment of gi(x). 
Depending on the interpretation on the discrepancy between experiment and the naive ex- 
pectation for (i.e., the Ellis- Jaffe sum rule [10]), two different sets of polarized parton 
distributions are often presented in the literature in the following way. First, one makes 
some parametrizations for spin-dependent parton densities based on some plausible (model) 
constraints. Then fitting these parametrizations to the data of gi(x) etc., one obtains (i) a 
best fit of Au(x), Ad(x) and As(x) at fixed Ql by assuming AG(x, Ql) = 0, or (ii) a best fit 
of AG(x) and the polarized valence distributions Au v (x), Ad v (x) with no sea polarization. 

However, most of the parton spin densities presented in the literature are problematic. 
First, model- independent QCD constraints on the valence spin densities Au v (x) and Ad v (x) 
at x — > 1 are not respected in many existing parametrizations. Second, most authors fail 
to employ a correct kernel A<7 lG (x), the hard cross section for photon-gluon scattering, to 
evaluate the gluonic contribution to the proton structure function g{{x). As we are going 
to stress in Sec. 2, whether or not gluons contribute to T\ is purely factorization dependent 
[11]. Once a factorization scheme is chosen, the "hard" kernel is completely fixed up to 
the factorization scale /if ac t- A determination of parton spin distributions using any other 
kernels, for instance the delta kernel, is certainly not trustworthy. 

In the present paper we shall give a critical analysis of the polarized parton distributions. 
We first give a brief overview in Section II on the role of the hard-gluonic contribution to the 
first moment of the polarized proton structure function. Based on the gauge-invariant and 
chiral-invariant factorization schemes, we then proceed to extract parton spin distributions 
from the g\ (x) data in Sections III and IV respectively. Sections V and VI contain discussions 
and conclusions. 

II. Framework 
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The sea-quark or anomalous gluonic interpretation for the violation of the Ellis- Jaffe sum 
rule depends on the factorization scheme defined for the quark spin density and the cross 
section for photon-gluon scattering. Much of the factorization scheme dependence and the 
related issues are already addressed by Bodwin and Qiu [11]. To set up the notation and 
the results necessary for our purposes, we will recapitulate the main points in Ref. [11]. 

The general expression of the proton structure function in the presence of QCD correc- 
tions to order a. is 
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where Af q depends on the regularization scheme chosen. Since the unpolarized parton 
distributions are usually parametrized and fitted to data in the MS scheme, it is natural to 
adopt the same regularization scheme for polarized parton distributions in which Af q (x) = 

fq( X ) ~ U 1 + X ) aIld ( See e -g-> [I 2 ]) 
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where the "+" distribution is given by 
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The first moment of f q {x) and Af q (x) is and —2 respectively. The parton spin densities 
in Eq.(l) are defined by Aq(x) = q^(x) +q^(x) — q^(x) —q^(x) and AG(x) = G^(x) — G^(x). 
It is known that a direct calculation of the polarized photon-gluon scattering box diagram 
indicates that Aa(x) has collinear and infrared singularities when m 2 = p 2 = 0, where m is 
the quark mass and p 2 is the four-momentum squared of the gluon. Depending on the choice 
of the soft cutoff, one obtains 



(i) m 2 = and p V [13] 
AoccmO) 




(4) 



(ii) m V and p 2 = [14] 



Aa AR (x) = (1 - 2x) f In ^ + In ^—^ 

\ m z x 



2(1- x) 



(5) 



1 It is known that /J Aa(x)dx = in the MS scheme [11]. However, the "hard" part of Aa(x) is dependent 
of the factorization scheme chosen, as elucidated below. For this reason, we shall discuss various soft cutoff 
schemes. 
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(iii) dimensional regularization [12] 

Aa R (x) = (1 - 2x) ( - + 7e + In + In ^—^ - 1 ] - 2(1 - x), (6) 

where /zms is a regulated scale in the minimal-subtraction scheme. For the first moment of 
Acr(x), it is easily seen that 

/ Acr CCM (x)Gfe = 1, [ Aa AR (x)dx = [ Aa R (x)dx = 0. (7) 
J o Jo Jo 

The result (7) can be understood as follows. Any term which is antisymmetric under x —>■ 1 — 
x, for instance terms proportional to const, x (1 — 2x), makes no contribution to Jq 1 Aa(x)dx, 
a consequence of chiral symmetry or helicity conservation, recalling that the gluon splitting 
function is of the form AP g c(x) = \{2x — 1). However, there is a chiral-symmetry-breaking 
term proportional to (1 — x) in the mass- regulator and dimensional regularization schemes, 
which compensates the hard contribution arising from the region k\ ~ Q 2 , where k± is the 
transverse momentum of the quark in the photon-gluon box diagram. 

Now, in order to consider hard-gluonic contributions to gf(x) (by "hard", we mean con- 
tributions with k\ > /if act ), one has to introduce a factorization scale /if act to subtract the 
unwanted soft contribution, i.e., the contribution arising from the distribution of quarks and 
antiquarks in a gluon: 

A^ rd (x, Q7/4J = Aa(x, Q 2 ) - Aa soi \x, tf act ). (8) 

In practice, one makes an approximate expression for the box diagrams that is valid for 
k\ « Q 2 and then introduces an ultraviolet cutoff on the integration variable k± to ensure 
that only the region k\ < /z 2 act contributes to the soft part [11]. The choice of the regulator 
specifies the factorization convention. When the ultraviolet cutoff is gauge invariant, it breaks 
chiral symmetry due to the presence of the axial anomaly and hence makes a contribution 
to Acr soft . Using the dimensional regulator for the ultraviolet cutoff it follows that [11] 

A^ccmW = (l-2x)(-i- 7E + ln^^ + ln-^-y-l^+2(l-x), 

= (l-2*)(---7E + ln^V (9) 
y e m z J 

Afff (i) = 0, 
for various soft cutoffs, and hence 

f 1 Aa^ u (x)dx = 1, / Ao%f£(x)dx= C Aa^\x)dx = 0. (10) 
Jo Jo Jo 



2 The last term —2(1 — x) in Eqs.(5) and (6) was neglected in the original work of Altarclli and Ross 
[14] and of Ratcliffe [12] respectively. It arises from chiral symmetry breaking due to the m 2 ^ cutoff in 
the mass- regulator scheme and the violation of the identity {7^, 75} = in the dimensional regularization 
scheme when e ^ 0. One may argue that this contribution is soft, for example, in the mass-regulator scheme 
if m 2 << /i 2 act and hence it does not contribute to "hard" Ac However, the cancellation of the ln(Q 2 /m 2 ) 
term, which depends logarithmically on the soft cutoff, from different x regions is not reliable because chiral 
symmetry may be broken at some hadronic scale. 
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In the mass-regulator and dimensional-regulator schemes, the original soft contributions in 
(5) and (6) are canceled by the contribution from chiral symmetry breaking introduced by 
the ultraviolet cutoff. Therefore, in the gauge-invariant factorization scheme 

Aa hard (x) = (1 - 2x) (in + In ^—^ - 1 ] - 2(1 - x), (11) 

where /i 2 act = 47r / u 2 4S exp(-7 E -l/e). It follows that $ Aff hard (x)rfi = 0. Note that Aa hard (x) 
is independent of the choice of soft and ultraviolet regulators. In this scheme, the quark spin 
has a gauge-invariant local operator definition: 

s^Aq = (plevMb), (12) 

where s M is the proton spin vector; it is Q 2 dependent because of the nonvanishing two- 
loop anomalous dimension associated with the flavor-singlet quark operator. The fact that 
gluons do not contribute to is in accordance with the OPE analysis in which only the 
quark operator contributes to T{ at the twist-2, spin-1 level [15]: 

jl = ^ (i - (^E^wsbV = \ (i - f ) + l Ad + l As ) >( 13 ) 

where Aq = Jq 1 Aq{x)dx. 

By contrast, it is also possible to choose a chiral-invariant but gauge- variant ultraviolet 
cutoff, so that [11] 

Q 2 , l-x \ „ , 2m 2 -p 2 x(l -2x) 

~2 i 2 2~T\ V + 111 1 - I 1 - ~2 , 2 2~R 

/j,f act + m z — p z x{± — x) x J /i£ act + m z — p z x{± 

and Jq 1 Aa ha,Td (x)dx = 1 for /i 2 >> p 2 , m 2 . In this chiral-invariant factorization scheme, the 
quark spin distributions in a gluon are obtained by a direct cutoff on the k± integration: 

Aq' G (x, /4J = f fact d 2 A; ± Ag G (x, fc ± ); (15) 

JO 

that is, all the quarks with k\ < /i 2 act in the gluon distribution are factored into the quark 
spin distribution. Contrary to the first scheme, Aq' = Jq 1 Aq'{x)dx cannot be written as a 
matrix element of a gauge-invariant local operator, but it is Q 2 independent as the gauge- 
variant ultraviolet cutoff in this scheme does not flip helicity; it is thus close and parallel 
to the naive intuition in the parton model that the quark helicity is not affected by gluon 
emissions. Replacing Aq(x) by Aq'(x) and Acr hard (x) by Acr hard (a;) in Eq.(l), we find 

jf 9l{x)dx = \ (l - ^) £ e 2 (Ag' - ^-AG). (16) 

3 Other main disparities between Aq and Aq' are as follows, (i) It is perhaps less known that [16] the 
spin-dependent Altarelli-Parisi evolution equations apply directly only to the gauge-invariant parton spin 
distributions. To evaluate the Q 2 evolution of Aq'(x) and Aq', one has to first apply Eq.(36) for example, 
(ii) In principle, Aq' and AG have a simple partonic definition: the former (latter) can be identified in 
one-jet (two-jet) events in polarized deep inelastic scattering [13]. 
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[1 - 2x) 



Consequently, Aq and Aq' are related by 

Aq = Aq' - ^-AG. (17) 

Finally we notice that it is also possible to choose an intermediate ultraviolet cutoff 
scheme which is neither gauge nor chiral invariant, so in general Aq = Aq' — X^AG for 
arbitrary A (A = and A = 1 corresponding to gauge- and chiral-invariant factorization 
schemes, respectively) [17]. It is clear that the issue of whether or not gluons contribute to 
Tf is purely a matter of the factorization scheme chosen in defining the quark spin density 
and the hard gluon-photon scattering cross section; a change of the factorization convention 
merely shifts the contribution of Aq(x) and Aa hard (x) in such a way that the physical proton- 
photon cross section remains unchanged. Though this controversy was resolved sometime 
ago by Bodwin and Qiu [11] (see also Manohar [18], Bass and Thomas [16]), it is consid- 
erably unfortunate that many of recent articles are still biased on the anomalous-gluonic 
interpretation of the g\ (x) data and that the work of Bodwin and Qiu is either overlooked 
or not widely recognized and well appreciated in the literature. 

III. Polarized parton distributions in the gauge-invariant factorization scheme 

This section is devoted to studying the spin-dependent valence and sea distributions 
based on the g\{x) data. We shall see that the positivity condition |As(x)| < s(x) due to 
the positivity of unpolarized parton distribution puts a very useful constraint on the shape of 
the polarized valence quark distributions. The presence of the gluon polarization will affect 
the shape of As(x), but not its first moment. 

To begin with, the combination of all EMC, SMC, E142 and E143 data for together 
with the SU(3) parameters [19] F + D = 1.2573 ± 0.0028 and 3F - D = 0.579 ± 0.026 yields 
[20] 

Au = 0.83 ±0.03, Ad= -0.43 ±0.03, As = -0.10 ± 0.03 , (18) 

and hence 

AS = Au + Ad + As = 0.31 ± 0.07 . (19) 

Decomposing Aq into its valence and sea components Aq = Aq v + Aq s , we shall follow 
Ref.[21] to assume that sea polarization is SU(3) invariant, i.e., Au s = Ad s = As. This 
assumption is justified since it leads to Au v + Ad v = 0.60 from Eq.(18), which is very close 
to the naive expectation that AS = 3F — D = 0.579 in the absence of sea polarization. The 
sea polarization is also found to be SU(3) symmetric within errors in the lattice calculation 
[8,9]. This is understandable since the disconnected insertions (for a definition of connected 
and disconnected insertions, see [8,9]), from which the sea-quark polarization originates, are 
presumably dominated by the triangle diagram and hence are independent of the light quark 
masses in the loop. (This effcet is absent in unpolarized distributions). Therefore, for SU(3) 
symmetric sea polarization, we obtain from (18) that 

Au v = 0.93 , Ad v = -0.33 . (20) 
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The Monte Carlo computation [8,9] shows that the magnitude of valence quark polarizations 
arising from the connected diagram is close to that given by (20). 



In terms of valence and sea spin distributions, Eq.(l) can be recast to the form 

(21) 



fj'lir.Q' 2 ) = ~ I -! ~Au v (y,Q 2 ) + ±Ad v (y,Q 2 ) + ^As(y,Q 2 





Recall that in the gauge-invariant factorization scheme, gluons contribute to g[(x), but not 
to rf. In general, both sea quarks and gluons contribute to the polarized structure function, 
but we will begin with the extreme case (i) As(x) ^ 0, AG(x) = 0. Since the unpolarized 
sea distribution is small at x > 0.2, the positivity constraint |As(x)| < s(x) implies that the 
data of gi(x) at x > 0.2 should be almost accounted for by Au v (x) and Ad v (x). Therefore, 
the shape of the spin-dependent valence quark densities is nicely restricted by the measured 
gi(x) at x > 0.2 together with the first-moment constraint (20) and the perturbative QCD 
requirement [22] that valence quarks at x — 1 remember the spin of the parent proton, i.e., 
Au v (x) I u v (x) , Ad v (x)/d v (x) — > 1 as x — > 1. In order to ensure the validity of the positivity 
condition |Ag„(x)| < q v {x), we choose the MRS (A') set [23] parametrized in the MS scheme 
at Q 2 = 4 GeV 2 as unpolarized valence parton distributions 

u v (x,Q 2 = 4GeV 2 ) = 2.26 x" a441 (l - x) 3 ' 96 (l - 0.54^+ 4.65x), 

d v (x,Q 2 = 4 GeV 2 ) = 0.279 x-°' 665 (l - x) 4 ' 46 (l + 6.80^ + 1.93x). (22) 

Accordingly, we must employ the same MS scheme for polarized parton distributions [see 
Eq.(l)] in order to apply the positivity constraint. For the spin- dependent valence distribu- 
tions we assume that they have the form 



Aq v (x) = x a (l- x) p {a + hyfx + cx + dx 1 ' 5 ), 



(23) 



with a and (3 given by Eq.(22). We find that an additional term proportional to x 1 ' 5 is 
needed in (23) in order to satisfy the above three constraints. 

For the data of gl(x), we will use the SMC [1] and EMC [6] results, both being measured 
at the mean value of Qq = 10 GeV 2 . Following the SMC analysis we have used the new -^(^O 
structure function measured by NMC [24], which has a better accuracy at low x, to update 
the EMC data (see Fig. 1). The best least \ 2 fit to g\(x) at x > 0.2 by (23) is found to be g 

Au v {x,Ql) = a;- a441 (l -a;) 3 - 96 (0.928 + 0.149^- 1.141a; + 11.612X 1 ' 5 ), 
Ad v {x,Ql) = ^-^(l-^^^-O.OSS-O^Sv^-S^eOx + S^Sx 1 - 5 ), (24) 



4 It was assumed in Ref.[21] that Au v (x) = a(x)u v (x), Ad v (x) = (5(x)d v {x) with a(x), (3(x) — > 1 as 
x — > 1 and a(x), [3{x) — > as x — > 0. However, the constraint at x = is not a consequence of QCD. In the 
present work we find that Au v (x)/u v (x) = 0.41 and Ad v (x)/d v (x) — —0.136 at x = 0. As a result, |Ag„(:r)| 
is usually larger than |As(x)| even at very small x (see Fig. 5 below). 
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at Qg = lOGeV 2 , which satisfies all aforementioned constraints. [] Since Ad v is negative 
while Ad v (x) is positive as x — > 1, it means that the sign of Ad v (x) flips somewhere at 
x = Xq [25]. We find that x = 0.496 in our case. 

It is evident from Fig. 1 that a negative sea polarization is required to explain the observed 
<7i(x) at small x. Assuming AG(x, Ql) = at this moment, we find from (21), (24) and the 
data of g{(x) that the polarized strange quark distribution is determined to be 

As{x, Ql) = -aT L17 (l - x) 963 (0.013Vz + 0.862x - 1.186x 15 ), (25) 

with As = -0.109 and x 2 /d.o.f. = 12.24/22, where uses of /i fact ~ 1 GeV and Q 2 = 
Qq = 10 GeV 2 have been made. It is easily seen from Fig. 2 that the positivity condition 
\As(x) / s(x)\ < 1 is respected. 

To illustrate the importance of having a least x 2 fit of g\{x) at x > 0.2 by Au v (x) and 
Ad v (x), let us consider another parametrization as an example | 

Au v (x) = 0.3588 x~ a54 (l -x) 3 - 64 (l + 18.36x), 

Ad v {x) = -0.1559 a;- a54 (l -x) 4 - 64 (l + 18.36a;), (26) 

with Au v = 0.93 and Ad v = —0.33 . It is evident that, contrary to Fig. 1, this parametriza- 
tion gives a reasonable eye-fit to the data (though x 2 /d.o.f. = 30/22) even at small x, 
as depicted in Fig. 3. One cannot tell if there is a truly discrepancy between theory and 
experiment unless the first moment of gi(x) is calculated and compared with data, i.e., 
(r?)theor y = 0.176 ± 0.006 versus (r?) expt = 0.142 ± 0.008 ± 0.011 [1]. Following the same 
procedure as before, we find that the sea polarization necessary to fit the data violates the 
positivity condition when x > 0.2 . This example gives a nice demonstration that an eye-fit 
to the data can be quite misleading. Therefore, we conclude that in order to satisfy the 
positivity constraint due to sea polarization, valence quark spin densities should fully accout 
for the observed gf(x) at x > 0.2. As a consequence, a deviation of theory from experiment 
for the polarized structure function should manifest at small x. 

It has been argued that a bound on As, namely \As\ < 0.052io;Q 23 , [27], can be de- 
rived based on the information of the behavior of s(x) measured in deep inelastic neutrino 
experiments and on the positivity constraint. However, this argument is quite controver- 
sial [28]. We note that since the strange quark distribution parametrized by MRS (A') [23] 
yields Jq xs(x)dx = 0.0182 for the strange sea momentum, it is consistent with the bound 
Jq xs(x)dx < 0.048 ± 0.022 extracted from the neutrino-nucleon experiment [29]. Hence, 
our As (a;) does satisfy all known constraints. More importantly, a sea polarization of order 
—0.11 in the polarized proton is confirmed by lattice calculations [8,9]. 

In a realistic case, it is very unlikely that AG(x, Q 2 ) vanishes at some scale Ql for all 
x. Even if AG(x, Qq) = at Q 2 = Ql, it can be radiatively generated at Q 2 > Q^. In the 
absence of any information on the shape and the magnitude of gluon polarization except for 

5 We have evoluted q v (x,Q 2 ) from Q 2 = 4 GeV 2 to 10 GeV 2 in order to compare with /S.q v {x,Q 2 t ). 
6 This parametrization is taken from Ref. [26] expect that we have made a different normalization in order 
to satisfy the first-moment constraint (20). 
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the restriction \AG(x)/G(x)\ < 1, we first take 

AG(x,Q 2 ) = 2.5A G (l-x) 7A \ (27) 

with Aq = 8.44 and AG = 2.5 , as an illustration. This parametrization is taken from the 
set A of gluon distribution in Ref.[25] but with a different normalization for our purpose. 
We see from Fig. 4 that the effect of polarized gluons is to suppress gi(x) at x < 0.01 and 
enhance gi(x) at 0.01 < x < 0.15 so that the net contribution to Y\ vanishes; that is, hard 
gluons contribute to g{{x) but not to in the gauge-invariant factorization scheme. Since a 
realistic polarized gluon distribution ought to have its first moment lie somewhere between 
and 2.5 , we take 

AG(x, Ql) = 0.199x" 117 (l - x) 5 - 33 (0.03 - 1.71y/x + 3.01a: + 43.5a; L5 ), (28) 

as determined below for case (iv). The first moment of this gluon spin density is AG = 0.5 . 
The presence of AG(x) will affect the shape of As(x) but not its first moment. Following 
the same extracting procedure as before for As(x), we find 

As(x, Ql) = -x~ l - 17 {l - x) 9m (0.0UVx + 0.865a; - 1.189a: 1 - 5 ), (29) 

with As = —0.11 and x 2 /d.o.i. = 11.75/22. The parametrizations (28) and (29) are regarded 
as the representative spin-dependent parton distributions for case (ii), as exhibited in Fig. 5. 



IV. Polarized parton distributions in the chiral-invariant factorization scheme 



As elaborated on in Sec. II, in the chiral-invariant factorization scheme the quark spin 
Aq' is Q 2 independent, and gluons contribute to the first moment of the polarized proton 
structure function. Since Aq' = Aq + f^AG and As v = 0, it is obvious that Aq' v = Aq v 
for SU(3) symmetric sea polarization. We shall follow Ref. [21] to assume that this is also 
true for their x dependence, i.e., Aq' v (x) = Aq v (x). Since AG(x) is also independent of the 
factorization chosen, we thus have 



9 P i(x,Q 2 ) 



x 




^Au v (y, Q 2 ) + ±Ad v (y, Q 2 ) + ^As'(y, Q 2 ) 



(30) 
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AG(y,Q 2 ) , 



with Aa hard being given by (14). We note that several different expressions for the kernel 
have been employed in the literature. For example, Aa(z) — 8(1 — z) was used by Altarelli 
and Stirling [30], Aa(z) = (1 — 2z) ln[(l — z)/z] by Ellis et al. [31] and by Ross and Roberts 
[32]. However, as we have stressed in Sec. II, a correct procedure of subtracting the soft 
contribution from Act (z) will yield a unique Aa haxA (z) up to the factorization scale //f ac t; 
which is independent of the choice of soft and ultraviolet regulators. 

We first discuss the extreme case, namely (iii) As'(a;) = and AG(x) ^ 0, which is just 
opposite to the other extreme case (i). It has been advocated that [33,13] a total absence 
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of sea polarization and an anomalous gluonic contribution might offer an attractive and 
plausible solution to the so-called "proton spin crisis" by accounting for the discrepancy 
between the Ellis- Jaffe sum rule and experiment for . It follows from (18) that AS' = 0.58 
with As' = 0, consistent with what expected from the relativistic quark model. To implement 
a large AE' and a vanishing As' demands a large gluon polarization: AG = — (27r/a s )As = 
2.5 at Q 2 = lOGeV 2 . The question then is: Can the data of g\{x) be explained solely by 
Au v (x), Ad v (x) and AG{x) without sea polarization? To examine this issue, we note that 
the gluon polarization is subject to the constraint 



J(x,Q 2 ) = j ^A^d (*,«-) AG(y,Q 2 ), (31) 
6tt Jx y \y ni ct J 

where 

J(x,Q 2 ) = g{{x,Q 2 



1 r 1 dy_ 
2Jy y 



^Au v (y,Q 2 ) + ^Ad v (y,Q 2 ) 



x 




(32) 



One may ask: Apart from the positivity constraints, can one treat Au v (x), Ad v (x) and 
AG(x) as free parameters and fit them to the measured g\{x)l The point is that when one 
works in the gauge-invariant factorization scheme, the shape of the polarized valence quark 
distributions, which is factorization scheme independent, is constrained by the positivity 
condition |As(a;)| < s(x), in particular in the region x > 0.2. Therefore, the l.h.s. of (31) is 
basically fixed by the data of g{{x) and the phenomenological Au v (x) and Ad v (x) given by 
(24). The polarized gluon distribution can be extracted from the Mellin transformation of 
(31) (for a detail of the procedure, see Ref.[21]). 



is 



The best least squares fit we found (see Fig. 6) for /i fact ~ 1 GeV and Q 2 = Ql = 10 GeV 2 

AG(x, Ql) = x- 1A7 (l - x) 5 - 33 (0.03 - 1.71<y/x + 3.01a; + 43.5x L5 ) (33) 

with x 2 /d.o.f. = 10.4/22 and AG = 2.51. There are two salient features with this spin- 
dependent gluon density: (i) AG(x) is negative at very small x, x < 0.025 . This is because 
the best x 2 fit to J(x) is positive at small x. We find that AG(x) = corresponds to 
a maximum xJ(x) occurred at x ~ 0.025 . Consequently, a negative behavior of AG(x) 
at very small x is natural, (ii) the positivity constraint |AG(rr)| < G(x) is violated at 
x > 0.15 (see Fig. 6). Hence, this AG(x) is physically unacceptable. However, we note 
that a fit to the g\{x) data with the polarized gluon distribution (27), which does respect 
the positivity condition, is equally acceptable with x 2 /d.o.f. = 14.13/22 (see the thick solid 
curve in Fig. 7). We thus conclude that it is still possible to reproduce the data of g\(x) 
with anomalous gluonic contributions (of course, the shape of the gluon spin distribution 
is basically arbitrary), yet a best least x 2 fit to data with x 2 /d.o.f. = 10.4/22 demands a 
polarized gluon distribution violating the positivity constraint. Needless to say, we have to 
await high-quality data in the future to pin down the issue. 
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There exist in the literature various parametrizations for polarized parton distribution 
functions fitted to the data within the framework of the chiral- invariant factorization scheme. 
However, most of them are not reliable or trustworthy owing to the incorrect use of the 
hard cross section Acr hard for photon-gluon scattering, among other things. For example, 
the predicted g\{x) using (26) for Aq v (x), (27) for AG(x) together with the delta kernel 
Aa hard (z) = 5(1 — z) fits the data very well with x 2 /d.o.f. = 11.9/22 (see the solid curve 
in Fig. 7). But the same set of parton spin distributions fails to fit the data at small x, 
x < 0.01, when the correct kernel (14) is employed (shown by the dotted curve in Fig. 7 
with x 2 /d.o.f. = 18.3/22). This is because gluon contributions at small x gain more weight 
via the convolution with the non-delta kernel. Recall that the same set of polarized valence 
quark distributions also leads to an unacceptable sea polarization when fitted to the data 
(see Sec. III). Therefore, we believe that our valence quark spin distributions parametrized 
by (24) are more sensible than any others. 

Since in a realistic case it is likely that sea polarization is nonvanishing and the value 
of AG is between and 2.5 , this leads to the more realistic case (iv) As'(x) ^ and 
AG(x) 7^ 0. If we assume that the shape of AG(x) remains the same as that of (33), the 
positivity condition of the gluon distribution requires that 

AG(x, Ql) = 0.199aT 117 (l - x) 5 - 33 (0.03 - 1.71y/x + 3.01a; + 43.5a; L5 ), (34) 

corresponding to AG = 0.5. Substituting this into Eq.(30) determines As'(x), which we 
find can be parametrized as 

As'(x, Ql) = -aT L17 (l - x) 9m (0M^/x + 0.69a; - 0.949a; L5 ) (35) 

with As' = -0.087 and x 2 /d.o.f. = 11.8/22 . 

V. Discussions 

In Sections III and IV we have considered four different cases for polarized parton dis- 
tributions in the gauge-invariant and chiral- invariant factorization schemes: (i) As(x) ^ 
0, AG(x) = 0, (ii) As(x) ^ 0, AG(x) ^ 0, (iii) As'(x) = 0, AG(x) ^ 0, and (iv) 
As'(a;) 7^ 0, AG(x) ^ 0. Since the value of AG ought to lie somewhere between and 
2.5 , it appears to us that case (ii) or case (iv) is more realistic. Note that cases (ii) and (iv) 
are not totally independent. This is because for a given AG(a;), which is factorization scheme 
independent, the sea quark spin distributions A<^(a;) and Aq s (x) are not independent and 
they are related via [21] (see also Eq.(57) of [16]) 

A&x) = Aq s (x) + * /' f 1 - -) AG(y), (36) 
7T Jx y \ yj 

derived from Eqs.(l), (11) and (14), where the assumption Aq' v (x) = Aq v (x) has been 
made (see Sec. IV). Clearly, its first moment is precisely Eq.(17), as it should be. We have 
explicitly checked that As(a;) of (29) and As'(a;) of (35) do satisfy the relation (36). The 
spin-dependent parton distributions in this work are presented in the MS scheme in the 
next-to-leading order of QCD (for a similar work, see [34]). 
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It is straightforward to compute the polarized structure functions g™ (x) of the neutron, 
gf(x) of the deuteron and their first moments and Tf respectively. The various polarized 
distributions satisfy the relation 

g{{x) + g^x) = 1 _l 5(jD 9i(z), (37) 
with ujd = 0.058 being the probability that the deutron is in a D state. We find 



= -0.053 , r? = 0.040 , at Q 2 = 10 GeV 2 , 



(38) 



while experimentally [1-6], 



■nn 
1 1 



-0.022 ±0.007 ±0.009, 
-0.037 ±0.008 ±0.011, 
{ -0.063 ±0.024 ±0.013, 



E142 at (Q 2 ) 
E143 at (Q 2 ) 
SMC at (Q 2 ) 



2 GeV 2 , 

3 GeV 2 , 
10 GeV 2 , 



(39) 



and 



■pa 
1 1 



0.034 ±0.009 ±0.006, SMC 
0.042 ±0.003 ±0.004, E143 



at (Q 2 ) = 10 GeV 2 , 
at (Q 2 ) = 3 GeV 2 . 



(40) 



Shown in Figs. 8 and 9 are the predicted polarized structure functions xg™(x) and xgf(x) 
respectively at Q 2 = 10 GeV 2 using the parton spin distributions in cases (i) and (ii). We see 
that although our predictions are consistent in gross with experiments, new measurements 
of g[ l (x) and gf(x) with refined accuracy are certainly needed. Note that the Q 2 dependence 
of polarized structure functions is not discussed here since our parton spin distributions are 
parametrized at Ql = 10 GeV 2 and it is known that only Q 2 > Ql evolution is governed by 
the Altarelli-Parisi equations. 



VI. Conclusions 



The fact that the size of the hard-gluonic contribution to = Jq g\dx is purely a matter 
of the factorization convention chosen in defining the quark spin distribution promotes us to 
consider four different possibilities of polarized parton distribution functions in two extreme 
factorization schemes: gauge-invariant and chiral-invariant ones. One cannot tell experimen- 
tally whether or not gluons contribute to Y\. We stressed that the hard cross section for 
photon-gluon scattering is unique up to the factorization scale //f act and is independent of 
the choice of the soft and ultraviolet regulators. 

Owing to the positivity constraints for sea and gluon polarizations, g\ix) at x > 0.2 
should receive almost all contributions from polarized valence quark distributions. This to- 
gether with the first moment and perturbative QCD constraints puts a very nice restriction 
on the shape of Au v (x) and Ad v (x). Eq.(24) is our best result for valence quark spin dis- 
tributions at average {Q 2 } = 10 GeV 2 . Working in the gauge-invariant factorization scheme, 
we have extracted the polarized sea distribution function from the EMC and SMC data of 
<7i(x) with the results (25) and (29) for cases (i) and (ii) respectively. All polarized parton 
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distributions in this work are presented in the next-to-leading order of QCD at the scale 
Q 2 = 10 GeV 2 . 

Based on the chiral invariant scheme and the aforementioned valence quark spin densities, 
we have found that it is possible to explain the measurements of g\{x) with anomalous 
gluonic contributions, yet a least x 2 fit to the data indicates that the best fitted gluon 
spin distribution violates the positivity condition |AG(a;)| < G(x). We have considered a 
more realistic set of parton spin distributions with a moderate value of AG = 0.5 and with 
a nonvanishing sea polarization. Many parametrizations of polarized parton distributions 
presented in the literature are not trustworthy due mainly to the use of an incorrect hard 
cross section for photon-gluon scattering. 

In principle, the choice of the set of Aq(x), AG(x), Aa hard (z) orofAg'(x), AG(x), Aa hard (z) 
to describe the polarized hadron structure function is a matter of convention. In fact, for 
a given AG(x), Aq'(x) and Aq(x) are related via Eq.(36). In practice, the gauge-invariant 
quantity Ag is probably more convenient and natural to use since it can be expressed as a 
nucleon matrix element of a local gauge-invariant operator. It is calculable in lattice QCD 
and, more importantly, its Q 2 evolution is directly governed by the polarized Altarelli-Parisi 
equations, which is not the case for Aq' and Aq'(x) [see the footnote after Eq.(15)]. 

Of course, inclusive polarized deep inelastic scattering experiments alone cannot reveal 
the magnitude and shape of the gluon spin distribution, and one has to await measurements 
of AG in independent processes in order to fully understand the proton spin structure. 
Nevertheless, there does exist a truly theoretical progress since the EMC measurement of g^, 
namely the lattice calculation of the proton matrix elements of the axial current [8,9]. The 
empirical SU(3) invariance observed by lattice QCD for the sea polarization manifested in 
the disconnected insertion strongly suggests that it is the axial anomaly which is responsible 
for the negative sea polarization and which explains the smallness of the quark spin content 
of the proton. 
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FIGURE CAPTIONS 



Fig. 1 The theoretical curve of xgf(x) fitted to the EMC and SMC data at x > 0.2 with 
the polarized valence quark distributions given by (24) and without sea and gluon 
polarizations. 

Fig. 2 The polarized strange quark distribution — As(x) fitted to the data of g\(x). Also 
shown is the unpolarized strange quark distribution evaluated at Q 2 = 10 GeV 2 using 
the MRS (A') parametrization [23]. 

Fig. 3 The predicted curve of xg\(x) arising from the spin-dependent valence quark spin 
distributions given by (26) without sea and gluon contributions. At first sight, it 
appears to give a reasonable eye-fit to the data. 

Fig. 4 Two theoretical curves for xg\(x). The solid line is the predicted curve for case (i) 
with x 2 /d.o.f. = 12.24/22, and the dotted curve with x 2 /d.o.f. = 14.95/22 is for case 
(i) plus the polarized gluon distribution given by (27). 

Fig. 5 Parton spin distributions for case (ii) parametrized at Q 2 = 10 GeV 2 . 

Fig. 6 The polarized gluon distribution extracted from a best least x 2 fit to the data of g{(x) 
by assuming As'(x) = 0. Also shown is the unpolarized gluon distribution evaluated 
at Q 2 = 10 GeV 2 using the MRS(A') parametrization [23]. 

Fig. 7 Three theoretical curves for xg[(x). With (27) for the polarized gluon distribution, 
the thick solid curve is calculated using (24) for valence quark spin distributions and 
(14) for the kernel, while the solid and dotted curves are based on (26) for Aq v (x) and 
the delta kernel Aa(z) = 5(1 — z) for the former curve and the kernel (14) for the 
latter. 

Fig. 8 The predicted polarized structure function gi(x) of the neutron at Q 2 = 10 GeV 2 
using the parton spin distributions in cases (i) and (ii). Also shown are the E142, E143 
and SMC data at the average Q 2 of each x bin. 

Fig. 9 Same as Fig. 8 except for the deuteron. The SMC data of gf(x) are evaluated at 
Q 2 = 10 GeV 2 , while E143 data at the average Q 2 of each x bin. 
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